Theory of composite-band Wannier states and order-N electronic-structure 
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From the order-N electronic-structure formulation, a Hamiltonian is derived, of which lowest eigen 
state is the generalized or composite-band Wannier state. This Hamiltonian maps the locality of 
the Wannier state to that of a virtual impurity state and to a perturbation from a bonding orbital. 
These theories are demonstrated in the diamond-structure solids, where the Wannier states are 
constructed by a practical order-N algorithm with the Hamiltonian. The results give a prototypical 
picture of the Wannier states in covalent-bonded systems. 

PACS number: 71.15.-m, 71.20.-b, 71.23.An 






a 



X3 

o 
> 

O 

On 
OS 
-I— > 



I 

O 

o 



X 



Recently, Wannier state (WS) has been re-focused as a 
foundation of the order-N methods, hnear-scahng meth- 
ods for large-scale electronic-structure calculations [0,|| . 
The original concept of WS's is defined within a single 
isolated band under the periodic boundary, but now is 
generalized to systems with composite bands and/or non- 
periodic boundaries, m^- This letter is devoted to the 
fundamental theories of WS's in composite-band systems. 

The generalized WS's can be defined, in insulators, as 
localized one-electron states that satisfy 



N 



H-ipk = y^gfcjV'j, 



(1) 



where N is the number of occupied states. Equation m) 
is derived from a variational procedure within a single 
Slater determinant. The parameters Sij are the Lagrange 
multipliers for the orthogonality constraints (V'ilV'j) = ^ij 
and satisfy Sji = {ipilHlipj). The above definition does 
not uniquely determine the wavefunctions. The resul- 
tant set of one-electron states {ipi} has a 'gauge' free- 
dom in the sense that any physical quantity is invari- 
ant under the unitary transforms with respect to the 
occupied states ipi ^ ijj'^ = '^j^iUijipj, where U is a 
unitary matrix. If this 'gauge' is fixed so as to diag- 
onalize the matrix Sij, we obtain the set of the eigen 



states {'0^'^''}' O'" ^^^ Bloch states in the periodic bound- 
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ary; H^,'^'^' = e^'^^Vr , where {e^^'} are the eigen 
energies. With WS's, the matrix e^ is not diagonal, 
but its trace gives the correct band-structure energy 

The diamond-structure solids, C, Si, Ge and a-Sn, 
are typical composite-band systems. For such materi- 
als, nearest-neighbor tight-binding (TB) Hamiltonians 
can be constructed within sp'^-hybridized orbitals, where 
the hopping along a bond is dominant. The correspond- 
ing hopping integral is the half of the difference between 
the energy level of an antibonding orbital (sa) and that 
of a bonding orbital (ffb); Aab = Sa,-^ £h- If all the other 
hoppings are ignored, the TB Hamiltonian is diagonal 



with respect to bonding and antibonding orbitals; 

N 

Ho = ^ (|bfe)eb(bfc| + |afe)£a(afc|) . 



(2) 



Here the k-th bonding and antibonding orbitals \hk), |afc) 
are defined by the pair of sp"^ orbitals on the A:-th 
bond. The WS's for Hq are just the bonding orbitals 
{|bfe)}fc=i,Ar, and this simple picture is the starting point 
of the present theory. In the present iterative calcula- 
tions of the WS's, the bonding orbitals are chosen as 
the initial states. Since all the bonds are symmetri- 
cally equivalent in the diamond structure, the resultant 
energy levels of the WS's {skk} have the unique value 

sws = (1/-^) X]i=i ^j 1 which is the weighted center of 
the valence band. 

Another important hopping in the diamond structure 
is the hopping within an atom, whose energy is one fourth 
of the energy difference between the atomic p- level (sp) 
and the s-level (Es); Aps = Sp — Eg- Within the present 
TB parameterizations, the electronic structures of the 
diamond-structure solids can be scaled with the unique 
parameter, am = Apg/Aab, called 'metallicity' and a sys- 
tem would be metallic, when am — > 1 |g|,^- For classifi- 
cation of group IV elements, some TB parameterizations 
were picked out and we obtained am — 0.44 for C 0, 
am = 0.75 [| for Si, and am = 0.77, for Ge §. In the 
present numerical demonstrations, we use the nearest- 
neighbor TB Hamiltonians H for Si whose parameters are 
from Ref. M. Here Aps is fixed to be 6.45 eV. Aab and all 
the other interatomic hoppings are functions of the bond 
length d. In the equilibrium case {d = do = 4.44a.u.), 
Aab = 8.25eV and am = 0.78. 

Equation ffl) is closely related to the localized-orbital 
order-N formulation g, where an energy functional 
^o(N) = Ef,(2<5,, - {iJjm){Hn\^j) = Tr[(2p- p2)r)] 
is iteratively minimized. Here f2 = H — -q, and p = 
X]fe=i lV'fe)('0/c|- P is the one-body density matrix and 
the energy parameter rj must be chosen to be sufficiently 
high {ri > En). The WS wavefunctions satisfy 



6E, 



■0(N) 



= {2n - pn - Qp)\^k) 



= 2(i/-iJocc)|^fe)-277(l-p)|^fe). (3) 

-ffocc is the Hamiltonian within the valence or occupied 
Hilbert space 



pH = Hp = Ho 



N 



)er\^r'\ 



(4) 



On the other hands, Eq.(P and the orthogonaUty 
constraint are rewritten as {H — -ffocc) IV"*:) = ^-nd 
(1 — p) iV'fc) = 0, which satisfy Eq.(||). In the present 
numerical calculations, the value rj = 5a. u. is chosen. 
From Eq. (|3|), we derive an eigen- value equation; 



^wslV'fe) =e/c/c|V'fc), 



where 



H^l = H- p^n - flpk 

N 

Pk^ p- |^fe)(V'fc| = E l^j)(^i 
j(#fe) 



(5) 

(6) 
(7) 



This eigen-value problem corresponds to the variational 
procedure of a specified WS (V'fc), while all the other 
WS's ({V'iljy/c) are fixed. If pk\i^k) = is satisfied, 
Eq.(p|) is reduced to Eq.(IT]) . A WS is not an eigen state 

of H, but an eigen state of ^-^s- ^^ '"^^^ ^^^ specified 
state (ipk) as 'central' WS. 

Figure |l| shows the density of state (DOS) of H and 
iJ^g in the Si case (d = do)- We see the follow- 
ing properties; (a) The 'central' WS {ipk) is the non- 
degenerate ground state of H^^g with the eigen value of 
£kk = ews- (b)The eigen states in the conduction band 
of H, (-0^ '^\ TV +1 <i< 27V), are also the eigen states 

r(fc) 



of -ffws with the same eigen energies 
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(c) All the other {N — 1) occupied WS's of H are not 

(k) 

eigen states of -ff^g- They form a high-energy band, lo- 
cated at e >2r] — ejy ~ 272 eV. The corresponding DOS 
profile is -Dhigh(£) ~ Dvai{2r] — e), where -Dvai(£) is the 
DOS profile of the valence band of H. This property rises 
from the relation 



(^,;|i?^^|V,) ^2rj- (^,;|i?occ|^,), I, J ^ k. 



(9) 



Due to the large energy shift of Dhigh(£), the Hilbert 
space of the other WS's is automatically 'excluded' in 
the variational freedom of the WS -ipk- 
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FIG. 1. The DOS of Hamiltonians H and H^^l for the Si 
case with 512 atoms in the periodic cell. The isolated level 
ews is broadened by a Gaussian with a width O.OleV and all 
the other levels are by that with a width of O.leV. The two 



parallel arrows in H indicate e\, and e 
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To construct WS's, we used an iterative order-N algo- 
rithm with the Hamiltonian -ffws • The periodic cell con- 
tains 4096 atoms and A^ — 8192 doubly-occupied WS's. 
For each WS -0^, the localization center was chosen at 
the center of the initial bonding orbital |bfe) and each WS 
was expanded into 614 sp^ orbitals with a spherical cut- 
off from the localization center. The resultant energy per 
WS (ews) has a deviation of about 0.006 eV (0.1%) from 
the correct value, where the correct value is obtained by a 
standard diagonalization method with the primitive cell 
and many k-points. The actual procedures are foUowings; 
(i) With proper initial states of WS's, the density matrix 
p and the Hamiltonians {ii^wg} are constructed, (ii) For 
each WS \4'k), the one-body energy (0fc|i?-^g|0fe) is min- 
(8) imized under the localization constraint on \^k) with the 

fixed Hamiltonian i/^s- ^™) "^^^ updated WS's {il^k}k 
are orthogonalized using the Lowdin symmetric orthog- 
onalization ||l^. Then the procedure goes back to (i), 
untill converges. Since the present TB Hamiltonians are 
upper-bounded, WS's can be also defined for the unoc- 
cupied or conduction band. The resultant WS's satisfy 
Eq.(|l]), where the N one-electron states should be those 
in the unoccupied band. Such conduction WS's can be 
formulated within the energy maximization procedure of 
Eo(^) , where the initial WS's is chosen as the antibond- 
ing orbitals {\a.k)}k=i.N and the energy parameter rj is 
chosen to be enough low. 

Figure || shows the norm distributions |Cfe0p of some 



WS's, where {4>} = {bfe,afc}. Figure || contains all the 
614 orbitals of the resultant WS's. The case (a) is the 
WS with d — O.Sdo, whose metallicity am — 0.47 might 
corresponds to the Carbon case. The case (b) is the equi- 
librium Si case (d = do)- For both cases, the norm dis- 
tributions on bonding orbitals are very small except the 
central one, because they are occupied mainly by the 
other WS's. The conduction WS's in the Si case is also 
plotted in Fig]|(c), which shows the similar decay prop- 
erty as in the valence WS (b), but the role of bonding 
and antibonding orbitals exchange with each other. The 
norm of the central bond is about 96 % in (a) or 94 % 
in (b) and (c). The summation of the norms upto the 
bondstep of n = 2 is about 99.8 or 99.7 % in all the 



cases. 
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FIG. 2. Norm distributions of WS's |Cfc^p, as a func- 
tion of the bondstep from the central bond. The closed cir- 
cles and the open squares denote the norms on bonding and 
anti-bonding orbitals, respectively, (a) the valence WS in 
the case with d — O.Sdo, (b) the valence WS in the Si case 
{d = do), (c) the conduction WS in the Si case. The crosses 
denote the values from the perturbation theory. In (a), the 
two values |C''^'''p and |C"^'p from the perturbation theory 
are almost identical. 



The most important issue for the practical order- 
N calculations is to reproduce physical quantities un- 
der localization constraints on WS's. With the exact 
WS's {V'fc}, tbe physical quantity of a one-body opera- 
tor X can be described as {X) = Ef (^i^^^l^lV'i"*^') = 
J2k {4'k\X\il!k)- In a practical order-N calculation with 
a localization constraint, the strict orthogonality con- 



straint is modified to an approximate one ((V'jIV'i) ~ ^ij) 
and the expression of a physical quantity is replaced by 
{X) = Ejj(2i5u - (V'j#»))(V'»l-^IV'i)7 which is mainly 
contributed by the diagonal elements {ipk\X\tpk)- If an 
operator X is short-ranged, like the present TB Hamil- 
tonian, its matrix element {ipk\X\^k) should be deter- 
mined dominantly by a 'central' region of WS, but little 
by the 'tail' or the asymptotic long-distance behavior. 
From this point of view, we analyse the present WS's 



through the matrix elements {ipk\X\ipk) of some opera- 
tors, not through the asymptotic long-distance behavior. 
Notice that the asymptotic long-distance behavior of WS 
is discussed in Refs. [Q. 

If we apply the above discussion to the operator 
(f — Tfe)^, where r^ is the center of the central bond 
(|bfc)), an effective spatial spread of WS can be defined 
as fws = {{''Pkli'f' — ■»"/c)^|V'fe))"^^^- The values were actu- 
ally calculated with the assumption that each sp'^ orbital 
is localized at the atom that the orbital belongs to. For 
a bonding orbital, this parameter is fb = d/2 from its 
definition. For the WS in the Si case {d = do), we ob- 
tained fws = 1.16rb. To see the effect of the boundary 
condition, we also calculated the system with 512 atoms 
without localization constraint and found fws = 1.19fb. 
Since the operator (f — Tfe)^ is not short-ranged, the value 
of fws rnight be sensitive to the boundary conditions. 

Another definition of the spatial spread of WS can be 



derived from the Hamiltonian H, 

rik) 



(k) 
WS- 



The Hamiltonian 



1 ~ £ws- A simplest case is the 
where the WS's are reduced to 



if^g has one localized eigen state i/'fe and the conduction 
band, and so H^^a, rnaps the WS formally to an impu- 
rity state, where the corresponding ionization energy is 
defined as Aws = ^X 
Hamiltonian in Eq.(p 
bonding orbitals and the corresponding gap Aws is to 
Aab. Using the uncertainty relation, a spatial spread is 
defined as ^b = 7i/-\/2'TioAab, where rric = la.u. Us- 
ing the value in the Si case (Aab =8.25eV), we obtain 
^b=0.29(io, which is consistent to the fact that the spread 
of a bonding orbital should be less than or about equal 
to the bond length d^. Such parameters can be also de- 
fined for WS's as Cws = ?i/\/2moAws = Cbv^Aab/Aws- 
For the WS in the Si case, we obtained Aws = 6.49eV 
and ^ws/fb — 1.13, which agrees, in the order, to the 
another definition of the spatial spread fws/^b — 1-16 
or 1.19. This agreement shows that the mapping theory 
to a virtual impurity state is consistent to the resultant 
WS's. The resultant values of the spatial spread lead us 
to the conclusion that the WS in the Si case is so local- 
ized that its spatial spread is in the same order as that 
of a bonding orbital. 

To see a limiting case with vanishing the bandgap, we 
also calculated an artificial case with Aps = 0, where 
the system is a direct-gap insulator and the bandgap 
A = SExx is located at T point ||T^ . We modified the pa- 
rameter Exx from the value in the Si case {E^x — 0.20eV) 
to an almost vanishing one {Exx = O.OOOSeV). This mod- 
ification was done by the tuning of Vppa and Vppjr , so as to 
keep Exy unchanged p2| . This modification changes the 
value of Aab as well. The resultant WS still shows a local- 
ized property with the spatial spread of fws/^b — 1.15 or 
1.20, where the latter value is from the calculation with 
512 atoms without localization constraint. From the gap 
parameters (Aab = 7.95eV,Aws = 4.23eV), we obtained 
Cws/Cb = 1-37. These resultant values of the spatial 



spread lead us to the same conclusions as those in the 
above Si case. 

The first-order perturbation of Eq.dq) can be con- 
structed using Hq in Eq.(l3) as the non-perturbative 
Hamiltonian; 



i^fc) = cW|bfe)+ J2 c^'^'^^ia, 



(10) 



Ji^^k) 



where C'-^-' ~ 1. The suffix v specifies the bond step and 
the inequivalent bond sites from the central bond |bfc). 
In the perturbation terms, bonding orbitals {|bj)}j^fc are 
'excluded', because these are the other WS in the non- 
perturbative terms and are in the high-energy band in 
Fig.|l|. Because the Hamiltonian H is a short-range op- 
erator, the perturbation series in Eq.([10|) contain only 
the 6 first-nearest-neighbor (FNN) antibonding orbitals 
(ja'-^-')) and the 18 second- nearest-neighbor (SNN) anti- 
bonding orbitals (ja'-^-')). For the FNN antibonding or- 
bitals, the perturbative coefficients are given ||] by 



C(o) 



x«|i/|fc;. 



a„ 



^ab 



8Aab 



(11) 



Here the factor 1/8 stems from the four atomic coordi- 
nation of the diamond structure, which is a three dimen- 
sional effect. On the other hand, the SNN anti-bonding 
orbitals are classified into two geometrically inequivalent 
bond sites; The 6 SNN bonds are parallel to the central 
bond. The other 12 SNN bonds exist, in a rough sense, 
within the plane perpendicular to the central bond. We 
denote the corresponding coefficients C'-^"' and C*-^^\ 
respectively, and propose the estimations of 



C(2A) (a(2^)|iJ|6fe) 



C(o) 



iab 



(12) 



where (2A) indicates (2 ||) or (2 _L). The first term is the 
first-order perturbation and its value is about -1-1/34 for 
C(2|t) oj. -1/27 for C(^^). Since this term is reduced to 
the ratio between two inter- atomic hoppings, its value 
is almost unchanged within the diamond-structure solids 
g. The second term in Eq.([l2|) is responsible for the 
successive hopping of the FNN hoppings, where C'-*'-' = 1 
is assumed. This term varies with am and is essential 
to the distinction between the cases in Fig.g (a) and 
(b). Notice that, in the realistic cases (am < 1), the 
second term (am/8)^ < 1/64 is smaller than the first 
term, though not negligible. The above coefficients can 
be determined by Eqs.(O), (|l2) and the normalization 
condition (|C(0)|2 + 6|C(^ + 6|C7(2ll)|2 + i2|C(2^)|2 = i). 

For the conduction WS's, the same perturbation theory 
can be constructed, where the role of bonding and an- 
tibonding orbitals exchange with each other. The re- 
suhant norm distributions jC^^^p, IC'^iDp and |C(2-L)|2 
are shown as crosses in Fig.0. The energy of the WS 
^ws = {'4'k\H\ipk) was also estimated from the pertur- 
bation results and the deviation from the correct value 



was about 0.06 eV in the Si case, which corresponds to 1 
% of the energy (ews) and to 10 % of the energy differ- 
ence from a bonding orbital {sh — £ws)- Here we can see 
that the present TB Hamiltonian is a short-range oper- 
ator and the value of its matrix element (ipklHlipk) can 
be well explained within a quite local area. 

In conclusion, the concept of composite-band WS's 
connects the picture of 'chemical bond' with the modern 
variational theory of electronic-structures. In this letter, 
we have shown how a WS is similar to and different from 
a bonding orbital within the diamond-structure solids, 
where the Hamiltonian i^ws plays a crucial role both for 
the construction and the analysis of WS's. These theo- 
ries are derived from the variational order-N formulation 
and so are applicable to other WS's in covalent-bonded 
systems and/or ab initio Hamiltonians. Results of these 
theories give microscopic pictures for practical order-N 
calculations of large-scale systems. 
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